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Abstract 

We have discussed the validity of the factorization approximation (FA) and nonextensivity- 
induced correlation, by using the multivariate (/-Gaussian probability distribution function (PDF) 
for iV-unit identical, independent nonextensive systems. The Tsallis entropy is shown to be ex- 
pressed by Sq = Sq N J A + ASg N ^ where q denotes the entropic index, S q N J A a contribution in the 
FA, and A.Sq a correction term. It is pointed out that the correction term of ASq N ^ is considerable 
for large \q— 1| and/or large N because the multivariate PDF cannot be expressed by the factorized 
form which is assumed in the FA. This implies that the pseudoadditivity of the Tsallis entropy, 
which is obtained with PDFs in the FA, does not hold although it is commonly postulated in the 
literatures. We have calculated correlations defined by C m = ((5xi 5xj) m ) q — ((5xi) m ) q ((5xj) m ) q 
for i ^ j, where 5x{ = Xi — {xi) q and {-) q stands for g-average over the escort PDF. It has been 
shown that C\ expresses the intrinsic correlation and that C m with m > 2 signifies correlation 
induced by nonextensivity whose physical origin is elucidated within the superstatistics. PDFs 
calculated for the classical ideal gas and harmonic oscillator are compared with the g-Gaussian 
PDF. A discussion on the g-product PDF is presented also. 

PACS numbers: 89.70. Cf, 05.70.-a, 05.10.Gg 
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I. INTRODUCTION 



In the last decade, much attention has been paid to the nonextensive statistics since 
Tsallis proposed the so-called Tsallis entropy The Tsallis entropy for iV-unit 

nonextensive systems is defined by 



s m = £-(i-J ] rfo lx) Yi x y (i) 

where q is the entropic index, ks the Boltzmann constant (ks = 1 hereafter), p q N \x) 
denotes the iV-variate probability distribution function (PDF), x = {xi} (i = 1 to N), and 
dx = YliLi d-Xi. The Tsallis entropy is a one-parameter generalization of the Boltzmann- 
Gibbs entropy, to which the Tsallis entropy reduces in the limit of q — > 1.0. The Tsallis 
entropy is nonextensive (non-additive), which is shown in the literatures as follows 
When the PDF for two independent subsystems A and B is factorized into those of A and 
B (x x e A, x 2 e B), 

pf{x u x 2 ) = p^(x 1 )p^(x 2 ), (2) 

Eq. (jU yields 

S®{A + B) = S®{A) + SU(B) + (1 - q)S^(A) S^(B), (3) 

which is referred to as the pseudoadditive relation. When the PDF for iV-unit independent 
subsystems is given as factorized form, 

N 



x) = HpV)(xi), (4) 



8=1 



we obtain the pseudoadditive Tsallis entropy Sq which is expressed by 



N 



In [1 + (1 - q)S< N } ] = ^ In [l + (1 - q)S^(z)] . (5) 

i=l 

It should be, however, noted that Eqs. ([3]) and ([SD are not correct in the strict sense because 
the bivariate PDF derived by the maximum-entropy method (MEM) cannot be expressed 
by Eq. (|2J) or (j3J), as will be shown shortly [Eq. (134"|) ]. Indeed, our calculation for identical, 
independent systems with the use of exact multivariate PDFs to Eq. ([I]) yields 

Sf = 2S« + (1 - q) [S^] 2 + ASf, (6) 
S™ = S$> A + ASW (7) 



where S q N J A denotes S q N ^ in Eq. (jSJ) evaluated by the PDF of Eq. (j4]) in the factorized 
approximation (FA), and ASq^ expresses a correction term [Eq. fpLtl) ]. Equations and 
show that does not satisfy the pseudoadditivity. 

The PDF is evaluated by the MEM for the Tsallis entropy with imposing some constraints. 
At the moment, there are four possible MEMs: (a) original method jl|, (b) un-normalized 
method 5/], (c) normalized method and (d) the optimal Lagrange multiplier (OLM) 
method fj. A comparison among the four MEMs is made in Ref. 4]. Although the four 
methods are equivalent in the sense that PDFs derived in them are easily transformed to 
each other [3] , obtained expressions for physical quantities are ostensibly different depending 
on the adopted MEM. 

Let us consider iV-unit independent systems whose hamiltonian is given by 



A' 



H = £ hi- (8) 



i=l 



PDFs for the hamiltonian in the Boltzmann-Gibbs statistics (q = 1.0) may be calculated 
with the use of either H or hi because they are given by 



N 



Pi 



{N) ( X ) oc Tr e~ pH = JJ p?{x t ), (9) 



i=l 



P P(xi) oc Tr,e^, (10) 

where (3 is the inverse of temperature, Tr denotes the full trace and Trj the partial trace 
over i. It is, however, not the case in the nonextensive statistics in which PDFs are given 
by 

p q N \x) oc Tr exp q (-PH), (11) 

JV 

^II^^) for 1.0, (12) 

i=i 

pWfe) oc Tri exp q (—phi). (13) 

Here exp q (x) expresses the g-exponential function defined by 

exp g (x) = [l + (l-g)x]V (1 - 9) , (14) 

with [x] + = max(x, 0), which reduces to exp g (a;) = e x for q — > 1.0. The inequality in Eq. 
( 1121) arises from the properties of the g-exponential function, 

exp q (x + y) ^ exp q (x)exp q (y) for q ^ 1.0. (15) 

3 



Then it has been controversial whether we should employ H or hi in calculating the PDF in 
the nonextensive statistics. In many applications of the nonextensive statistics, one usually 
calculates Pi(xj) with the use of hi in Eq. (1T51) . explicitly or implicitly employing the FA, 
because an exact evaluation of Pn{ x ) i n Eq. ffTTl) is generally difficult. This issue of the 
degree of freedom N on the PDF has been discussed in Refs. js, 9), [hJ. A calculation for 
TV-unit harmonic oscillator shows that the partition function obtained in the FA is quite 



different from that obtained by the exact PDF 
IIIC. 

In our previous papers 



111 ], related discussion being given in Sec. 



12 



131 ]. we discussed the effect of spatial correlation on the Tsallis 
entropy and the generalized Fisher information in nonextensive systems. We obtained the 
multivariate g-Gaussian PDF with the OLM-MEM [6|, which correctly includes correlation. 
It is the purpose of the present paper to discuss the issue mentioned above, by using the exact 
multivariate g-Gaussian PDF derived in Ref. [12]]. One of the advantages of a use of the 
g-Gaussian PDF is that it is free from an ambiguity in defining the physical temperature 
in conformity with the zeroth law of thermodynamics in the nonextensive statistics [4j]. 
From calculations of the Tsallis entropy and correlation in iV-unit independent systems, 
we will show the importance of effects which are not taken into account in the FA. The 
nonextensivity-induced correlation has been discussed for classical ideal gas Q, Q and 
harmonic oscillator [151 ] . 

The superstatistics is one of alternative approaches to the nonextensive statistics besides 



n n 



jy a superposition over the fluctuating 
13]- [2^]. The superstatistics has been 



the MEM jl7|, Il8|, Il9] (for a recent review, see [2£j]). In the superstatistics, it is assumed 
that locally the equilibrium state of a given system is described by the Boltzmann-Gibbs 
statistics and its global properties may be expressed 
intensive parameter (i.e., the inverse temperature) 
adopted in many kinds of subjects such as hydrodynamic turbulence, cosmic ray and solar 
flares 20j. The physical origin of the nonextensivity-induced correlation may be elucidated 
within the superstatistics. 

The paper is organized as follows. In Sec. II, multivariate PDFs for correlated nonex- 
tensive systems derived by the OLM-MEM jo] are briefly discussed [12I . [ljj. By using the 
multivariate PDF, we calculate the Tsallis entropy and correlations. Some model calcu- 
lations of the q- and iV-dependent Tsallis entropy and correlations are presented. In Sec. 
Ill, we discuss the physical origin of the nonextensivity-induced correlation, calculating the 



4 



PDF within the superstatistics [l?j, Q]- PDFs of the one- dimensional classical ideal gas 
and harmonic oscillators derived with the use of the OLM-MEM [6j are compared with the 
g-Gaussian PDF. The PDF expressed by the g-product [21] is also discussed. Sec. IV is 
devoted to our conclusion. 



II. 4-GAUSSIAN PDF 



A. OLM-MEM 



We consider iV-unit nonextensive systems whose PDF, pj^Hx), is derived with the use of 
the OLM-MEM [(3] for the Tsallis entropy given by Eq. ([1]) [l, 2]. We impose four constraints 



given by (for details, see Appendix B of Ref. 

1 = jp q N \x)dx, 



i 

i=l 
1 N 

i=l 

. N N 



(16) 
(17) 

(18) 

(19) 



i=l j=l(&) 



Here /i, a 2 and s express the mean, variance, and degree of intrinsic correlation, respectively, 
and (-) q denotes the g-average over the escort PDF, 



Pf)(*) 



Pq N \x) 



Cq 



with 



I [pi N \x)] q dx. 



(20) 



(21) 



Evaluations of the g-average with the use of the exact approach 
Appendix A. 



22 



231 ] are discussed in 



The OLM-MEM with the constraints given by Eqs. f[T6|) - f[T9|) leads to the PDF given by 



12] 



Pq N \x) 



7 {N) eX P<? 



<f>(x) 



(22) 



where 



N N 



$(a;) 



a 
01 



^ y^jao Sij + ai(l - - - ft), 

i=i i=i 

[l + (iV-2) 3 ] 

(l- S )[l + (iV-l) S ]' 
s 



;i- s )[i + (iv-i) s ] 



(23) 

(24) 
(25) 



27TI/, 



Wo 



9-1 

rf^vra 2 )^ 2 



iV/2 r /_i j\n 



, qr— 1 2 



.(AO 



2nu, 



1 N / 2 



l-q 



r(i^+i+f; 



for g > 1, 
for g = 1, 
for q < 1, 



(26) 



{( l_ s) ^l [ l + (iV _l )s]} l/2 ) 

(JV + 2) - iVg 



(27) 
(28) 



B(x, y) and r(z) denoting the beta and gamma functions, respectively. Hereafter we assume 
that the entropic index q takes a value, 

2 



< q < 1 



iV' 



(29) 



because p q (x) given by Eq. fl22|) has the probability properties with i/g > for q < 
1 + 2/N and because the Tsallis entropy is stable for q > 24]. 
In the absence of the intrinsic correlation (s = 0) for which 



(30) 



Eq. (|22D reduces to 



,(iV) ex Pg 



1 



N 



Wo J i=l 



(31) 



On the other hand, the PDF in the FA is given by 



N 



p q N FAW = n p?m> 

1=1 

1 N 



exp. 



(Xi - flf 



(32) 
(33) 
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In Eqs. (j3~Ij) and fl33|) . Z\ is given by Eq. fl26l) with n = 1.0. From a comparison 
between Eqs. (13T|) and (133|) . it is evident that 

P? ^) ^ (34) 

except for q = 1.0 or N = 1. 



B. Tsallis entropy 

Substituting the PDF given by Eqs. (j22])-(r58D to Eq. ([I]), we first calculate the Tsallis 



entropy which is given by 



1 - d N) 



q-l 

with 



= W* (35) 



r (JV) _ ,,(iv) 



[Zf)] 1 - 9 . (36) 



The s-dependence of the Tsallis entropy was previously discussed (see Fig. 1 of Ref. |12|). 
With increasing s, the Tsallis entropy is decreased as given by 

5f )( a ) = 5f >(0) - ^ ~ 1)C ^ s 2 for |«| « 2/y/N(N-l) . (37) 

Now we pay our attention to identical, independent systems with s = 0, for which we 
obtain 

Si N) = Siy A + AS^\ (38) 



with 



S { q % = (39) 



ASM = (^ T )[(4 1) ) Ar -cf)], (40) 

= (^t) ([^C^) 1 -^ " • (41) 

Here S^ N J A denotes the Tsallis entropy calculated in the FA, and ASg N ^ signifies a correction 
term. Equation (jSJ) leads to 

&A = £ W ^l^ 1 -<r\S?y, (42) 

= ^a) + ^-^-^) (g ay + ... (48) 



In particular for iV = 2, Eq. (138|) becomes 



with 



1 _ ^ 

S? = - i ^ r = S^ A + AS?\ (44) 



S$ A = = + (1 - g)(5«) 2 , (45) 



ASf = (J-lKcWf-a, (46) 



l 

g ^ T )(W 1) (^ 1) ) 1 -f-^ ) (^ 2) ) 1 - 9 )- (47) 

We will show some model calculations of S g N ^ and ASq N \ The q dependences of 
\og 10 (Sq N ^ / N) for various N values are shown in Fig. 1(a) where the inset shows an en- 
larged plot for q > 1.0. With decreasing q from unity, S q N ^ is logarithmically increased. 
The N dependence of ASq / Sq N ^ is shown in Fig. 1(b) where the inset again shows an 
enlarged plot for q > 1.0. We note that AS^/S^ is 0.049, 0.0 and 0.025 for q = 0.0, 1.0 
and 2.0, respectively. When N is more increased, AS q N ^ / Sq NS> becomes more considerable. 
For example, its value becomes 0.61 and 0.91 for N = 10 and 20, respectively, at q = 0.0. 

The N dependence of \og 10 (Sq NS> / N) is shown in Fig. 2(a). We note that with increasing 
N, Sq N ^ is increased (decreased) for q < 1.0 (q > 1.0). The N dependence of ASq / Sq N ^ is 
shown in Fig. 2(b) where results for q = 0.9, 1.1 and 1.5 are multiplied by a factor of ten. 
It is realized that for q = 0.2 and 0.5, ASq N ^ / Sq^ — > 1.0 as iV — > 100 where S^ A almost 
completely underestimates the Tsallis entropy. 

Figures 1(b) and 2(b) clearly show that ASq N ^ / Sq is positive and becomes appreciable 
for large values of N and/or large \q — 1|, in particular for q < 1. 



C. Correlation 

Next we calculate correlations for i ^ j defined by 

C m = ((Sxi 5x 3 ) m ) q - ((SXtDq ((SXiDq, (48) 

where 5xi — Xi — //. With the use of the PDF given by Eqs. (!22~ l) -( l28l) . we obtain the first- 
and second-order correlations given by (for details, see Appendix A) 

d = a 2 s, (49) 
C 2 = C 2s + C 2n , (50) 

8 




FIG. 1: (Color online) The q dependence of (a) Sf'/N and (b) ASf } /Sf } for iV = 2 (chain 
curves), 4 (dotted curves), 10 (dashed curves) and 20 (solid curves): insets show enlarged plots for 
q > 1.0: the ordinates in (a) are in the logarithmic scale. 
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FIG. 2: (Color online) (a) log 10 [Sf } /N] and (b) ASf'/Sf 1 as a function of log 10 N for q = 0.2 
(solid curves), 0.5 (dashed curves), 0.9 (chain curves), 1.1 (squares) and 1.5 (circles): results for 
q = 0.9, 1.1 and 1.5 in (b) are multiplied by a factor often. Note that, for q = 1.0, Sq /N = 1.4189 
and ASq N ^ / Sq N ^ = independently of N. 



10 



where 



C, 



2s 



a 



2/i 



2[(iV + 2) -Nq] a A s 2 
[(N + 4) - (N + 2)q] ' 

for s = 0, 
oo for q -> 1 + 2/(iV + 2), 
2s 2 for g = 1.0 or AT — > oo, 

2(g-l)a 4 
[(JV + 4) - (N + 2)qY 

for g -> 1.0, 
oo forg-^l + 2/(iV + 2), 



— jr for N — > oo. 



(51) 



(52) 



(53) 



(54) 



We note that Ci and C*2 S arise from intrinsic correlation s for < q < 1 + 2/(iV + 2). In 
contrast, C*2 n expresses correlation induced by nonextensivity, which vanishes for q — 1.0 
and which approaches — 2/N as iV — ► oo. In particular for iV = 2, and C2„ are given by 

2(2 - q)a 4 s 2 



C*2n 



(3 - 2g) 
(q - 1) a 4 



(3 - 2q) ' 

On the contrary, the factorized PDF given by Eq. ( |33l) yields 



(55) 
(56) 



/~<FA _ r~iFA _ f~iF A _ n 



(57) 



By using Eq. MM, we obtain correlation of C m for arbitrary m with s = 0, 



A, 

/ 1 , 



2^ a 2 
(9-1) 



T{q/{q-l)~N/2-m) _ ( T{q/ {q-l)-N /2-m/2) 
r( 3 /( 3 -l)-JV/2) I r(g/(g-l)-2V/2) 



(1-9) 



r( g /(i-g)+jy/2+i) 

r(g/(l~g)+A r /2+m+l) 



r(g/(l-g)+jV/2+l) 
r(g/(l-g)+JV/2+m/2+l) 



for g > 1.0, 



for g < 1.0, 



(58) 



where 



A, 



r(l/2+m/2) 

r(i/2) 



for even m, 



for odd m. 



(59) 
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FIG. 3: (Color online) (a) The q dependence of C% s (the dashed curve), C 2n (the dashed curve) 
and C 2 {=C2s + C2n, the solid curve) for s = 0.5 and N = 2. (b) The q dependence of C2 with 
s = for N = 2 (the solid curve), 4 (the dotted curve), 10 (the dotted curve) and 20 (the chain 
curve) . 

The nonextensivity yields higher-order correlation of C m than m > 2 for q ^ 1.0 in inde- 
pendent nonextensive systems where s — C\ — 0. 

Figure 3(a) shows the q dependence of C 2s , C 2n and C 2 (= C 2s + C 2n ) for s = 0.5 and 
N = 2. We note that C 2s > 0, while C 2n < for q < 1.0 and C 2n > for q > 1.0. 
The g dependences of C 2 for s = 0.0 with iV = 2, 4, 10 and 20 are plotted in Fig.3(b), 
where C 2 vanishes for q = 1.0, and it is -0.33, -0.25, -0.14 and -0.083 for iV = 2, 4, 10 
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FIG. 4: (Color online) (a) The q dependence of C m with N = 2 and (b) the N dependence of C m 
with q = 0.5 (s = 0): m = 2 (the solid curve), 4 (the dashed curve) and 6 (the chain curve). 

and 20, respectively, at q — 0.0. Figure 4(a) and 4(b) show the q- and iV-dependent C m , 
respectively, for m = 2, 4 and 6 with s = 0.0. Magnitudes of C m are significantly increased 
with increasing m. 
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III. DISCUSSION 



A. PDF in the superstatistics 



The physical origin of the correlation induced by nonextensivity is easily understood in 
the superstatistics. We consider the iV-unit Langevin model subjected to additive noise 
given by [12j 



dx 



—1 = -\ Xi + V2D &(t) +1 for i = 1 to JV, 



(60) 



where A denotes the relaxation rate, £j(£) the white Gaussian noise with the intensity D, 
and / an external input. The PDF of n^ N '(x) for the system is given by 



N 



n- (i) 



8=1 



where the univariate PDF of Tt^'ixi) obeys the Fokker-Planck equation, 







[(Xxi - I)TX {l \x u t)} + D^^{x h t) 



2 



dt dxi 
The stationary PDF of n^\xi) is given by 

1 



dxi 



(1). 



V2 



: exp 



{Xj - (J,)' 
2a 2 



with 



(61) 



(62) 



(63) 



(64) 



After the concept in the superstatistics [17], [l8|, [l9|, |20j , we assume that a model parameter 



fi = I/X, a 2 = D/X. 



°f P (= JU -QX fluctuates, and that its distribution is expressed by the ^-distribution with 



rank n 



mm 



n / 2-1 e -n(3/ 20o 



(65) 



r(n/2) \2{3 Qj 

where T(x) is the gamma function. Average and variance of (3 are given by (/3)s = fio and 
({/3 2 )p-Po)/(3o = 2/n, respectively Taking the average of 7r^(ac) over /(/?), we obtain the 
stationary PDF given by 



13| 



AN) 



exp 



0o 



N 



2u, 



(TV) 



(66) 
(67) 



i=i 
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with 



q 



(N) "I N l 2 N 



1 + 



(g-i)A) 

2 



i=l 



2'g 



(iV + n) ' 



(68) 
(69) 



where v q is given by Eq. ( |28|) . In the limit of n — >■ oo (q — > 1.0) where /(/?) — >■ — /? ), 
the PDF reduces to the multivariate Gaussian distribution given by 

N 



(3 



H exp 



8=1 



(70) 



which agrees with Eqs. (IBTj) and (1551) for /3 = A/D = 1/cr 2 . 

We note that the PDF given by Eq. (167|) is equivalent to that given by Eq. (13T|) derived 
by the MEM when we read /3q = 1/a 2 . The nonextensivility-induced correlation arises 



from the common fluctuating field of (3, because pq(x) ^ FL^^fa) f° r 1 7^ 1-0 despite 



7T 



X 



(1): 



B. PDF for classical ideal gas 



It is worthwhile to discuss the PDF of one- dimensional ideal gas, whose hamiltonian is 
given by 

N N o 



8=1 8=1 

m and pi standing for the mass and momentum, respectively, of ideal gas. By employing 
the OLM-MEM [6], we obtain the PDF given by (for details, see Appendix B) 



r(N) 



exp„ 



v. 



(N) 



8=1 



2m 



(72) 



with 



(N) 1 N/2 p/^ N\ 

2irmu^ L \q-l 2j 



(?-l)/3 

N/2 



27r m 

/3 



for 1 < q < 3, 
for g = 1, 



(73) 



2irm v, 



(N) 



^ J for g < 1. 
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The internal energy is given by U = N/2/3 [Eq. (1B7I) ] independently of q, which yields the 
Dulong-Petit specific heat as in the Boltzmann-Gibbs statistics. The PDF of nonextensive 



ideal gas was originally discussed in Ref . |25j with the use of the normalized MEM . Later 
it was re-examined by the OLM-MEM 6] [26]. The PDF given by Eq. (|72|) is equivalent 
with the g-Gaussian PDF given by Eq. (1311) . if we read m/(3 — > a 2 . 



PDF for classical harmonic oscillator 



Next we consider the one-dimensional iV-unit harmonic oscillators whose hamiltonian is 
given by 



N 



N 



,2^.2 



(74) 

i=l i=l v 7 

m, uj, Xi and Pi expressing the mass, oscillator frequency, position and momentum, respec- 
tively. With the use of the OLM-MEM [6j, the PDF is given by (for details, see Appendix 
C) 



P ( q N) (p,x) 



r(N) 



exp. 




(75) 



with 



n (2N) - 
2-7T V q 


N 


r (^) 








N N 






»P) 






2nu q 2N ^ 


N 




(l-q)uj(3_ 




r(x^+i+iv) 



for 1 < g < 3, 
for g = 1, 
for g < 1, 



(76) 



where 



(2JV) 



(JV + 1) - iVg [Eq. 



The internal energy is given by U — N//3 [Eq. 



IT|) ]. which is the same as that in [15j and in the Boltzmann-Gibbs statistics. The PDF of 
the harmonic oscillator was calculated by using un-normalized [5] and normalized 2] MEMs 



in Ref. jllf. where U has rather complicated q and N dependences (see Eqs. (12) and (22) 
of Ref. |llj). 

PDFs for classical ideal gas [Eqs. (1721) ] and harmonic oscillator [Eq. (j75p ] have the same 
structure as the g-Gaussian PDF given by Eq. (|3"T1) . and they are expected to have the 



same properties as the g-Gaussian PDF. Actually correlation definec 
(i 7^ j) is shown to be induced by nonextensivity in ideal gas fl4 . 
oscillator 



3 QU I 



; (hihj)q- (hi) q (hj) q 
15l . 1 1 61 ] and harmonic 
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D. g-product PDF 

From functional forms of PDFs in Eq. (1B3I) or (1C3I) . we expect that p q (x) may be 
expressed by 

pf\x) = p«K) ® q p[ 1 \x 2 ) ® 9 • • ® q p q x \x N ), (77) 



where the g-product (<g) g ) is defined by [21] 



= [x 1 "^!/ 1 - 9 - l] 1/(1 - g) . (78) 

Equation (177)) . however, does not hold when we take into account the precise form of PDFs 
including their normalization factors. For example, for a univariate PDF given by 

pW(x) = -^[1 - (1 - g)^] 1 ^), (79) 

Xg 

with 

X® = Tn (80 ) 
Eqs. Q7_7J) and (J7JI]) yield the g-product PDF for N = 2, 

pf , (^i)® 9 pf ) (^) = , y(J L [2 - ff) 1 - 9 - (1 - + fe 2 )] V(1 ' g) ■ (81) 

(X q >) l ~1 

Unfortunately, Eq. (I8T1) is not in agreement with the bivariate PDF given by 

pf\ Xl , x 2 ) = -JL[1 - (1 - + M] VM , (82) 

pPM^pPM, (83) 

with 

= Tr 12 [1 - (1 - q)P{h x + h 2 )] 1 ^\ (84) 

It is easy to see that the relation: p q \xi, x 2 ) = p q (xi) ® q p q l \x 2 ) holds only for X^ = 
X ( q 2) = 1 in Eqs. (jSIl) and p^j). 
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IV. CONCLUDING REMARKS 



It has been postulated that the Tsallis entropy Sq satisfies the pseudoadditivity given 
by Eq. (j3D or (jSJ), and that the factorized PDF: p q N \x) = YliPv ( x i) leads to the Tsallis 
entropy expressed by Eq. ([T]) [3]. The pseudoadditivity is a basis of the Tsallis entropy 



and detailed discussions on its pseudoadditivity have been made 



27, 28 



29j . We should 



note, however, that this is not self-consistent because the PDF derived by the MEM for S, 



(N) 



leads to p q (x) ^ YiiPQ } { x i) which contradicts with the postulated, factorized PDF. The 
pseudoadditivity cannot be a basis of the Tsallis entropy 

It has been also controversial whether the statistical independence is expressed by (a) 



p q N \x) = YliPq 1 ' \ x i) or (b) H = hi in nonextensive systems 
the condition (a) expresses the statistical independence of iV-unit subsystems, we obtain 



If we assume that 



S, 



(at) _ gvy^ which satisfies the pseudoadditivity. However, the PDF derived by the MEM 
for Sf° yields p q N \x) ^ J] 

i p q 1 \x i ), which is inconsistent with the assumption. This means 
that the statistical independence cannot be expressed by neither the product nor g-product 
PDFs. 

Our calculations with the use of the multivariate g-Gaussian PDF have shown that 



(i) the Tsallis entropy is given by S, 



(AO 



S 



(AT) 
q,FA 



+ ASq N) > S, 



g,FA 



where the correction term 



of ASq NS> is significant for large \q — 1| and/or large N, 

(ii) the Tsallis entropy does not satisfy pseudoadditivity, and 

(iii) nonextensivity-induced correlation is realized in higher-order correlations C m for m > 2 
while C\ expresses the intrinsic correlation. 

The items (i)-(iii) are expected to hold also for classical ideal gas and harmonic oscillator. 
The item (ii) is against the common wisdom [30]. The nonextensivity-induced correlation 
in the item (iii) is elucidated as arising from common fluctuating field introduced in the 



superstatistics 



3,Q, 



19| . It has been shown that the FA is not a good approximating 



method in c 
in Refs. 31 



assical nonextensive systems, just as in quantum ones as recently pointed out 



321 ] . We should be careful in adopting the FA, although it has been widely 



employed in many applications of the classical and quantum nonextensive statistics [33]. 



18 



Acknowledgments 



This work is partly supported by a Grant-in-Aid for Scientific Research from the Japanese 
Ministry of Education, Culture, Sports, Science and Technology. 



APPENDIX: A. EVALUATIONS OF gr- AVERAGES 

We briefly discuss evaluations of g-averages given by 
#(a) = / [1 - (l-q)a$(x)} 1/{1 - q) dx, 



(Al) 



Qf )(a) = (Q(x)) 



1 



9 AN) 7 (N) 



Q(x) [1 - (1 - q)a $(x)] q/( - 1 - q) dx, (A2) 



by using the exact expressions for the gamma function 

■y poo 



22 



23 



3l|: 



y 



r(«) Jo 



u 



s - l e - yu du 



for s > 0, 



y s = — r(s + 1) / (-r a_1 )e-^ (it for s > 0, 
2?r y c 



(A3) 
(A4) 



where a = l/(2i/J JV) tr 2 ), $ (cc) is given by Eq. (123]) . Q(cc) denotes an arbitrary function 
of a;, C the Hankel path in the complex plane, and Eq. (I36I) being employed. We obtain 

2d, bd, lai] 



VTA Jo°° u^~ l e~ u z[ N) [(q - l)au\ du for q > 1.0, 

/ c (-*) -T ^ -1 e _ ^i (iV) [-(l - * for g < 1.0, 



r 



(A5) 



r 



/^^^"^-"^[(g - l)cra] <2?°[(g - 1)H for g > 1.0, 
/ (-t)~^-W° [-(!-?)«*] 



A + 1 



x Qf°[-(1 -g)at] tit 



for g < 1.0, 
(A6) 
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where 

Z [ N \a) = I e- a *Wdx, (A7) 



Q[ N \a) = I Q{x)e-"*Wdx. (A8) 

Z\ '{a) J 

Thus we may evaluate the q average of Q(x) from its average over the Gaussian PDF. 
For example, by using the relations for N = 2, 

<(^-^) 2 >i = ^ (A9) 
{(xi - n){xj ~ fj))i = 7T- for % ^ j (A10) 



2a 
_+ 

4a 2 



1 + 2s 2 

{(xi -n) 2 (xj - jj) 2 )i = — -=— fori^j, (All) 



and employing Eqs. (IA5I) and (IA6I) . we obtain 

((^-//) 2 > 9 = (A12) 
((a* - n)(xj - fi)) q = a 2 s for i ^ j, (A13) 

which yield Eqs. (jI5])-([5T]) and ([53]). 



APPENDIX: B. PDF FOR IDEAL GAS 



In order to obtain the PDF of pq N ^ (p) for classical ideal gas with the OLM-MEM [6( , we 
impose the constraints given by 



Jp { q N) (p)d P , 



The OLM-MEM |6| yields 



U ( 2m ) 

i=l N ' < 



(Bl) 
(B2) 



pf(p) oc 



iV 



1 V(l-9) 



i-(i-g)/? E 



i=l 



2m 
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(B3) 



where the Lagrange multiplier f3 expresses the inverse of temperature. Rewriting Eq. (1B3P 

as 



pM(p) oc 



1=1 



Pi 

2m 



1/(1-9) 



(B4) 
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with 



a 



g 

l + (l-g)/3Ef 



we obtain [/ in terms of a as given by 



17 



N 



(N + 2- Nq)a 

From Eqs. (1B5|) and (1B6|) . C/ and a are self-consistently determined as 

N 



U 



a 



2/3'- 



2(3 



N + 2-Nq 

With the use of Eqs. (jBijh flBTl) and (|B8]), the PDF is given by Eqs. ([72} and ([73 



(B5) 

(B6) 

(B7) 
(B8) 



APPENDIX: C. PDF FOR HARMONIC OSCILLATORS 



We obtain the PDF for classical harmonic oscillators with the OLM-MEM imposing 
the constraints given by 



1 



J J P ( q N) (Pi x ) dpdx, 



N 



The OLM-MEM [6| leads to 



i=l 



2^,2 



^\2m 2 



P ( q N) (P,x) oc 



which is rewritten as 



N 



i=l 



2m 2 



i V(i-3) 



P ( q N) (P,x) OC 



N / 1 2 2 \ 



V(l-9) 



with 



a 



/3 



We obtain U in terms of a as given by 

U = 



l + (l-q)/3U' 



N 

(N +1 — Nq)a>' 



(CI) 
(C2) 



(C3) 



(C4) 



(C5) 



(C6) 



21 



From Eqs. (1C5I) and (1C6j) . U and a are self-consistently determined as 

N 

U = j, (C7) 

» = FTiW < C8 > 

By using Eqs. fjCijl . flCTfjl and (EE}, we finally obtain the PDF given by Eqs. (1751) and (1751). 
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